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SUMMARY 


A common analysis procedure at the Goddard Space Flight Center is to analyze a 
structure to predict responses that will occur during a sinusoidal or random 
vibration test in which a predetermined motion is input by the shaker at the 
shaker/structure interface. Often the loads produced in these tests are a 
controlling factor in the design of the structure. The response analyses are 
generally performed in two stages: a modal analysis to determine the undamped 

' eigenvalues and eigenvectors of the structure "cantilevered" at the interface 
where the structure will attach to the shaker, followed by a subsequent steady 
state frequency response analysis using the modal information. The output of the 
modal analysis, while giving information that is useful in comparing with the 
■ eigenvalues and eigenvectors obtained in the test, does not give a sufficiently 
clear indication of which eigenvectors will be important contributors in the 
subsequent frequency response analysis. The reason for this is that the 
vibration test loading is really a distributed inertial loading on the v/hole 
structure known in terms of the input acceleration as opposed to some known force 
at a given location in other forced response analyses. 


By using the eigenvector data obtained in the normal mode analyses, a modal 
participation factor can be calculated which gives a clear indication of the 
relative importance of each calculated mode in terms of its response to any base 
acceleration input. This factor, like the generalized mass or generalized stiff- 
ness already calculated in NASTRAN, is a property o^ the structure. In addition, 
however, it is also a property of the base excitation. 

This paper presents a technique for calculating the modal participation factors 
for base excitation problems using a DMAP alter to the NASTRAN real eigenvalue 
analysis Rigid Format. The DMAP program automates the generation of the 
"seismic" mass to add to the degrees of freedom representing the shaker input 
directions and calculates the modal participation factors. These are shown in 
the paper to be a good measure of the maximum acceleration expected at any point 
on the structure when the subsequent frequency response analysis is run. 
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Structures that are analyzed to predict what their responses will bt in a 
vibration test quite often have one feature in common; namely that they are 
rigidly attached to a vibration shaker at a number of interface degrees of 
freedom. Such situation, therefore, can be analyzed by multi-point constraining 
all of the interface degrees of freedom to any convenient single interface grid 
point with specification of the shaker motion at this grid point. The shaker 
motion, or base excitation, as it were, is then specified for a nonredundant set 
of degrees of freedom (DOF) that can easily be identified on a NASTRAN SUPORT 
bulk data card. In the following derivation of the base excitation modal 
participation factors, it is assumed that the structure to be analyzed is to have 
a base motion at a nonredundant set of DOF and that the interface OOF have been 
constrained to this set of nonredundant DOF. In addition, the derivation is 
restricted to systems with real eigenvectors. 


DERIVATION OF MODAL PARTICIPATION FACTORS 


To derive an equation for the modal participation factors due to base excitation 
consider a general structural model, as shown in Figure 1. Presumably this 
structure would be rigidly connected to the shaker at several points, denoted in 
this example as 1 through 4. This situation can be simulated by creating some 
central point, r, and connecting 1 through 4 to r via rigid elements. Then r is 
a set of nonredundant DOF at which the shaker motion can be specified. The 
equations of motion of this structure, using the displacement set notation of 
NASTRAN, and ignoring deunping for the moment, are * 



Where contains the loads necessary to drive the 
motions 'imparted by the shaker. 


r 


( 1 ) 


DOF by the specified base 


Following elimination of multi-point constraints (via MPC or rigid element cards) , 
single point constraints (SPC's) and DOF omitted via Guyan reduction (OMIT's),^ 
in equation (1) is reduced to and the equations are 


where 








( 2 ) 


( 3 ) 


♦Bold face letters will be used to designate matrices, including row and column 
vectors . 
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The lAr- DOF are identified on SUPORT bulk data carda and 
DOF in Partitioning equation (2): 



are all remaining 



where are the forces exerted by the shaker on the structure and result in the 
base motions IA ^ . The forces ^ ^ are not kno«m , a priori } rather the motions 
(accelerations, velocities or displacements) at the r DOF are known. In this 
context the forces in Pj, are really forces of constraint which are necessary to 
produce the desired base motion. Notice in equac .on (4) that there are no applied 
loads for the Jl DOF. The situation being simulated here is one in which the 
only excitation comes eibout due to the shaker driving the structure at the »' DOF. 


The DOF Call be considered 


to be made up to two terms: 


where 





(5) 



displacements of the x DOF when the DOF are fixed 
rigid body displacements of the DOF due to the DOF 


Thus would be the displacement DOF for a structure that was "cantilevered" or 

fixed at the '' DOF. The displacements, at the a DOF due to the »' DOF not 

being fixed can be determined simply from geometric considerations and is 



0 U 


v\ 


(6) 


where matrix Oir data block DM generated by module RBMG3 in NASTRAN whenever 
SUPORT bulk data cards defining the existence of r DOF is present. The columns of 
are rigid body displacements at the ^ DOF due to a unit motion at one of the r 
DOF. They are rigid body displacements due to the fact that the F DOF are a 
nonredundant set. Substituting (5) and (6) into (4) the two matrix equations for 
the unknowns \A'^ (since Q^^is known) are 

V ' rtu )U, (7, 


M U' 1 

u 1 




( 8 ) 



0 







i 

f 

S 








9 


0 


where use of the fact that 


and thus 
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) 


( 9 ) 


has been made. This is discussed in the NASTRAN Theoretical Manual, Reference 1. 
Equation (7) is the well known form for base acceleration input where the problem 
is solved for the motions as if the base were fixed but with an equivalent 

inertia loading on each of the DOF that is the rigid tody inertia force due to 

the base moving. Equation (8) merely solves for the forces required to drive the 

base at the specified motion. A better form for P*. results if equations (9) and 

(7) are used to eliminate the stiffness term in (8) . The result is 


P. { ‘ Mj u, , n, u 


( 10 ) 




where 


M 0' W 0. * D’ M V 0 I H 


( 11 ) 


is the rigid body mass matrix of the structure relative to the SUPORT DOF and is 
data block MR generated by NASTRAN in module RBMG4 when SUPORT bulk data cards are 
used. Equation (7) can be diagonalized with a modal transformation using the 
"cantilevered" modes of the structure (that is, the modes of the structure when 
constrained at the r DOF) . Define 


A ' ‘ . 

= J eigenvector of structure when constrained at r DOF. 

I = matrix whose columns arc eigenvectors^ ^ 


win *•*- 

Then the modal transformation is 








( 12 ) 
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where j is a vector of modal DOF s', (1 =1, 2,... h ) with h being the number of 

modes used in the dynamic analysis. Substituting (12) into (7) and 


premultiplying by X"*" 

the result is 







•n.4. ' 

- ■ 


(13) 

where 

^ * 




- 

I 

Xin 

(14) 

and 

n 

T ^ tv ■ 

r 



9 K * 

Xv'i ’'u ‘ 

L. 

(15) 


are bcth diagonal matrices with m containing the generalized masses (denoted 
as (lf) j ) and generalized stiffnesses (denoted as fc j ) of the 

cantilevered modes. 


With the addition of damping, each individual modal equation in equation (13) is, 
.after dividing through by the generalized mass; 

f ’ ,T 

b * M„) U, (16) 


The damping factor, , is specified through use of the TABDMP bulk data card. 

Define 


.T 


1 


4' (m 0 4 M, ) 


(17) 


' The row matrix Xy. contains as many terms as there are DOF identified on the 
SUPORT card, that is, those that will be given a known motion from the shaker. 

; Each term in is identified as being a modal participation factor since the 

solution to equation (16) for some specified base motion at one DOF will be 

proportional to the corresponding term in * 

' 1 

As an example of this consider the situation when the base motion is a steady 
* state sinusoidal input at some resonant frequency tv • with an acceleration 
(complex) magnitude . In this case: ^ 


t 


u 






11 7 


s 
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The solution to equation (16) for the modal acceleration is 

^ I 


where 


); 


X ' L 




u 


1 


( 18 ) 


and 




I 


is the modal amplification in mode \ . Equation (18) indicates the modal 
acceleration is, of course, lagging the base acceleration input by 90 degrees. 

The significance of the modal participation factors in is evident from 

equation (18). At resonance, the modal acceleration magnitude is equal to the 
product of the modal amplification, the participation factor and the magnitude of 
the input shaker acceleration. For systems with relatively small damping, the 
magnitude of the accelerations at grid points throughout the structure can be 
approximated by a single mode response (the resonating mode) so that, for example, 
using only one mode in equation (12) 



(19) 


The modal participation factors, together with the eigenvectors calculated in a 
real eigenvalue analysis can therefore be used to obtain estimates of what the 
magnitude of the grid point accelerations will be in a subsequent forced base 
response analysis. The modal participation factors can be easily calculated 
during the real eigenvalue analysis run (in Rigid Format 3, for example). A 
matrix of all participation factors for all modes is: 




W i ( M D ' ^ y, ) 


( 20 ) 


That is, each row of is an ^ (the participation factors for one mode). 

Each column of p gives the participation factors for all modes for one DOF of 
base acceleration. For example, column 1 of contains the participation 

factors for base acceleration in the first DOF entered on the SUPORT card. 
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The modal participation factors can be used to easily calculate what is commonly 
referred to as modal mass. Consider equation (10) for the base loads (the shaker 
forces required to produce )• Using the modal transformation in equation (12) 
the loads in equation (10) can be written as 


I' 




where 


* ab • 

Using equation (18) , the complex amplitude of ^ is P where 

p; = pjc'“' 


Using equation (17) ; 


= Q- (22) 

where 

is the mode modal mass matrix. It is a square matrix of size r and from 
equations (21) and (22) it is seen that its use is in determining the j mode 
contribution to the base forces (required shaker forces) . For systems, with small 
damping the forces at a resonance u)j would be approximately so that 

equation (22) would give an approximation to the total forces. It is 
interesting to note that if the sum of all modal mass matrices is denoted as 


0 


z V. 


0 
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and if all of the eigenvectors are used (h-H.) then; 

If there are no masses at any of the V- set of D.O.F then 

'hv =• M (26) 

Yt V-y 

That is, the sum of the modal masses is the rigid body mass matrix (data block MR 
from module RBMG4) when all masses are at the r set of DOF. Thus, modal masses 
can be calculated in a real eigenvalue analysis and can be used to estimate how 
many modes will have to be used in the subsequent base excitation dynamic analysis 
to make sure a high percentage of the complete system mass is included (since 
modal truncation, in effect, truncates mass). This would be useful in transient 
analysis where the single mode response at resonance is not a valid approximation. 
The modal mass matrix can easily be determined from the modal participa- 
tion factors From equations (17) and (23) it is seen that 

'('if. 4! (27) 


DMAP ALTER DESCRIPTION 


The DMAP Alter package for calculating the modal participation factors and modal 
masses is described in detail in Appendix A. The calculations are done, of 
course, in the real eigenvalue analysis NASTRAN Rigid Format and in particular in 
either COSMIC Rigid Format 3 or MSC Rigid Formats 3 or 25. The description in 
Appendix A is an alter package for MSC Rigid Format 25 or COSMIC Rigid Format 3. 
However, calculation of the modal masses in the COSMIC version is not convenient 
at this time due to a limitation in the COSMIC version of module DIAGONAL. 

As recommended in the MSC documentation. Reference 2, base motion problems are 
solved by putting all DOF wh h are to have a specified motion on SUPORT cards 
and also connecting fictitious large masses to these DOF. If the large masses 
are very large in comparison to the complete structure mass then, when a force is 
applied to these DOF, their acceleration will be determined to be very nearly the 
force divided by the large mass. Making the force equal to the large mass tines 
the desired acceleration at these DOF will then result in the NASTRAN solution for 
these accelerations to be the desired base accelerations. 

Thus, the first portion of the DMAP Alter automates the inclusion of the large 
(seismic) massgs at the DOF identified on SUPORT bulk data cards. The large mass 
is taken as 10 times the diagonal of the rigid body mass matrix (data 
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block MR calculated in module RBMG4) . The bulk data deck for the real eigenvalue 
analysis with the DMAP Alter program should include SUPORT cards for the 
nonredundant DOF at which base motion is desired. As mentioned previously, the 
actual base DOF that are connected to the shaker are to be rigidly constrained to 
some nonredundant set of DOF that are then included on SUPORT cards, (and not 
single point constrained) to identify them as the base motion DOF. 

The second part of the DMAP Alter calculates the matrix which contains the 

modal participation factors. These are printed out as matrix data block FNR. 

The third -part of the DMAP alter calculates the diagonal terms of the modal mass 
matrix These can only be calculated in the MSC version since they require 

1 feature in the DIAGONAL module (namely raising all elements of a matrix to a 
power) not available in ,the COSMIC DIAGONAL module. The matrix which contains 
the diagonal terms of is data block MJRR and row ] of this matrix contains 

the diagonal terms for mode ] . This matrix is also printed out in terms of 
percent of the corresponding mass in the rigid body mass matrix as matrix data 
block MJRRP. Addition of each column should give a value of 100 if all modes are 
used and if there are no masses (other than the fictitious large mass) at the 
set of DOF. 


SAMPLE PROBLEM 


Figure 2 shows a sample problem using the DMAP Alter for a beam which is to be 
given a base motion at grid point 11. The aluminum beam is 2.54m (100 inches) 
long with a mass of 9.1kg (beam weight^i| 201b) ^ Section properties for the beam 
are: . 0013m^ (2. Oin'^) area and 8.33X10 m (2.0in ) moment of inertia. The JOF which 
are to be given base motion should be the higher numbered DOF in the problem to 
avoid numerical problems with the large fictitious masses added to these DOF. 

The fact that the seunple problem initially has a nonredundant interface with the 
shaker (at grid point 11) is of no consequence. As mentioned earlier, if the 
interface where the structure is to be connected to the shaker is redundant, the 
interface points can be rigidly connected to some one point that is then called 
the base motion point. 

Appendix B contains selected pages from the NASTRAN run for this sample problem 
including the bulk data deck listing. The input data is seen to resemble any 
real eigenvalue run for the beam cantilevered modes except that some of the DOF 
at the cantilever point (grid point 11) are identified on a SUPORT bulk data card 
instead of constrained. For the sample problem these were taken to be DOF 1, 3 
and 5 signifying that these will be the DOF that would be given a base motion in 
subsequent dynamic analyses. The remaining DOF at grid point 11 are constrained. 

As mentioned previously, the DMAP Alter automates the inclusion of the large 
fictitious masses at the DOF identified in the T set. The values of these masses 
are contained in matrix data block MAAB output in the DgAP Alter and included in 
Appendix B. As seen from MAAB, the large masses are 10 times the rigid body 
masses output from the grid point weight generator (taking into account the latter 
is in weight units and MAAB is in mass units) . 
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The eigenvalues for the sample problem are included in Appendix B. The first 
three are rigid body inodes due to the three DOF on the SUFORT card. The remaining 
modes are those that would have been found if the SUPORT DOF had been constrained 
instead. This is due to the fact that large masses are connected to the SUPORT 
DOF in order that the elastic modes will be the same as if these DOF were con- 
strained. Several of the modes are shown on Figure 3. 

The modal participation factors are shown in Appendix B as the matrix data block 
FNR. These are also shown in Table 1. Each column in the t 2 d)le is a column of 
and there are 3 columns) one each for the SUPORT DOF. Each row in the tcdjle ''' 
shows the terms in (which are the rows ) which are the participation 

factors for the mode corresponding to that row of the table. 

The second column of Table 1 gives the participation factors for base motion in 
DOF 3 at grid point 11; namely lateral shaker motion. From Table 1 it is seen 
that mode 1, with a frequency of lO.Shz, has a participation factor of 1.56 while 
mode 5 does not contribute at all. This is due to the fact that mode 5, as seen 
from Figure 3 is an uncoupled longitudinal mode of the beam and could have been so 
identified by visual inspection of the mode shapes for this simple beam problem. 
For more complicated structures the participation factors are more helpful in 
identifying significant modes for various base excitations than the visual 
inspection of the eigenvector data. 

The use of the participation factors in estimating grid point accelerations that 
would be found in a steady state sinusoidal shaker simulation is demonstrated on 
Figure 4. The figure shows comparison oJ data from an actual frequency response 
run using Rigid Format 8 along with the estimated accelerations based on the 
modal participation factors in table 1. 

The solid line is thu acceleration at the be 2 un tip as a function of the shaker 
frequency, , with nn input acceleration in the Z lateral direction of l.Sg's at 
grid point 11. The simple calculation using the participation factors and a 
single mode response with equation (19) gives acceleration magnitudes indicated 
by the data points on Figure 4. The good agreement shows the value of knowing 
the participation factors when the modal analysis is run rather than waiting to 
determine the frequency response in a later run. The data point for the 490. 6hz 
mode on Figure 4 is not a mode in resonance since the 490.6 mode is a longi- 
tudinal mode and the base response in the example is lateral. 

The output in Appendix B also shows the diagonals of the modal mass matrices 
as output matrix ^-JRR. Expressed as a percent of he corresponding mass term in 
the complete structure rigid body mass matrix, the result is output matrix 

MJRRP. As mentioned previously, ac there are no masses at the DOF (other than 
the fictitious large mass put there for use calculating the participation 
factors) and if all modes of the structure are used, the columns of MJRRP should 
add to 100. Table 2 shows the terms from MJRRP. As with the modal participation 
factors these also clearly indicate which modes are significant in terms of base 
mocion in certain directions. For general base excitation problems, the data 
from MJRRP is useful in identifying significant modes to be included in the 
analysis. The columns of Table 2 do not add up to 100 percent for base motion 
DOF 1 or 3 since the set did have mass. The beams' mass was distributed to all 
eleven grid points so that the end grid points (1 and 11) each have 5 percent of 
the total. 
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CONCLUSIONS 


Knowledge of the modal participation factors for base motion can be useful in 
determining the significant modes to be used in forced base motion problems. 
Obtaining the participation factors when the modal analysis is run, allows 
estimates to be made of the subsequent strrrture responses to th' ,se motion 
prior to running the base motion dynamic analysis. The particiiation factors use 
for steady state base motion f nquency response analyses was demonstrated; 
however they are equally as valid fcr transient analysis. 
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Table 1 

Modal Participation Factors 
for Bea..i Problem 


Modal Participation Factors From 
for Base Acceleration at G.P.ll In: 


Mode 

No. 

Frequency 

(Hz) 

DOF-1 

DOF -3 

DOF -5 

4 

10.94 

0. 

1.5*=C^ 

-113.590 

5 

67.82 

0. 

-.8446 

17.800 

6 

188.0 

0. 

.4736 

-6.124 

7 

364.5 

0. 

-.3337 

2.923 

8 

490.6 

1.2706 

0. 

0. 

9 

595.9 

0. 

.2161 

• 1.590 

10 

878.2 

0. 

- ,1593 

.9801 

11 

1202. 

0. 

.137-' 

-.7370 

12 

1460. 

-.4165 

0. 

0. 

13 

1545. 

0. 

.1154 

-.5638 

14 

1861. 

0. 

-.0806 

.3669 

15 

2086. 

0. 

-.0453 

.1987 

16 

2393. 

-.2414 

0. 

0. 

17 

3267. 

-.1632 

0. 

0. 

18 

4061. 

.1171 

0. 

0. 

19 

4755. 

-.0854 

0, 

0. 

20 

5332. 

.0813 

0. 

0. 

21 

5777. 

-.0414 

0. 

0. 

2? 

6081. 

.0240 

0. 

0. 

23 

6234. 

-.0079 

0. 

0. 


r 
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Table 2 



Modal Masses 
for Beam Problem 




Diagonals of Modal Mass Matrices* 
for Base Acceleration at G.P.ll In; 


ftode 

No. 

DOF-1 

DOF- 3 

DOF- 5 

4 

0. 

61.073 

97.030 

5 

0. 

18.854 

2.4995 

6 

0. 

6.4685 

. 3228 

7 

0. 

3.3013 

.0856 

8 

80.724 

0. 

C. 

9 

0. 

1.9882 

.0321 

10 

0 

1.3149 

.0149 

11 

0. 

.9087 

.0078 

12 

8.6749 

0. 

0. 

13 

0. 

.6166 

.0044 

14 

0. 

.3585 

.0022 

15 

0. 

.1171 

.0007 

16 

2.9142 

0. 

0. 

17 

1.3315 

0. 

0. 

18 

.6854 

0. 

0. 

19 

.3647 

0. 

0. 

20 

.1878 

0. 

0. 

21 

.0858 

0. 

0. 

22 

.0288 

0. 

0. 

23 

.0031 

0. 

0. 

Totals 

95.0002 

95.0008 

100.000 


‘Expressed as percent of total rigid body mass of structure. 
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General Base Excitation Prcfclan 


z 

4 
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xll0 98 76 54 32 1 

( r) 

Figure 2 

Sanple Beam Problem 
With Base Excitation At Grid Point 11 
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Jt3de 1 

Rigid Body x Translation 


Mode 2 

Rigid Body z Translation 


X 


Mode 3 Mode 4 

Rigid Body y Rotation 1st Bending, 10.94hz 



Mode 5 Mode 6 



Mode 7 



ftode 8 

1st Lcxigitudinal, 490. 6hz 


Figure 3 
Beam Mode Shapes 
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100. 



Frequency (hz) 

Figure 4 

Acceleration at G.P. 1 of Beam due to 
1.5g Sinusoidal Base Acceleration at G.P. 11 
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DMAP Alter for Determining Base Excitation ftodal Participation Factors 
(Real Eigenvalue Analysis Rigid Formats) 


1 

ALTER a 


2 

PARAM ^ 

//C,N,NOPA,N,T=-l $ 

3 

PARAMR 

//C,N,DIV/V,N,GRAVO/C,N,1.0A»Y,vm«ASS*“1.0 $ 

4 

PARAMR 

//c,n,complex/v,n,outr/v,n,gravo/c,n,o.o/v»n,grav $ 

5 

DIAGONAL 

MR/MRD/C , N , SQUARE $ 

6 

ADD 

MRD,AJRD/V,N,GRAV $ 

7 

ADD 

MRD,/MBASE/C,N, (l.E8,0.) $ 

8 

VEC 

USET/CP/C,N,A/C,N,L/C,N,R $ 

9 

MERGE, 

,,,MBASE,CP,/MAAB/C,N,-1 $ 

10 

MATGPR 

GPL,USET,SIL,MAAB//C,N,A/C,N,A $ 

11 

ADD 

MAAB,MAA/MAAT $ 

12 

EQUIV 

MAAT,MAA/T $ 

13 

ADD 

MR,MBASE/NENMR $ 

14 

EQUIV 

NEWMR.MR/T $ 

15 

ALTER a 


16 

PARTN ^ 

PHIA, , CP/PHIL, PHIR, ,/C,N,l/C,N,2/C,N,2/C,N,2 $ 

17 

GKAM, 

, PHIL, MI, LAMA, ,,, ,CASECC/MGNN, ,KGNN,PHILN/C,N, -1/C, N,-0/ 
C,Y,LFREQ=. 000001/C, Y,HFREQ=1.E20/C,N, -1/C, N, -1/C, N,-l/ 

v,n,nocupA»n,fmode $ 

18 

ADD 

MGNN,/WGNN/V,N,GRAV $ 

19 

MPYAD 

MLL,DM,MLR/ML3/C,N,0/C,N,1/C,N,0/C,N,2 $ 

20 

ADD 

ML3,/WL3A,N,GRAV $ 

21 

MPYAD 

PHILN,WL3,/FNR0/C,N,1/C,N,1/C,N,0/C,N,2 $ 

22 

SOLVE 

WGNN,FNR0/FNR/C,N,1 $ 

23 

MATPRN 

FNR,,,,// $ 

24 

DIAGCM^AL 

FNR/MJRRO/C, N, WHOLE/C, N, 2.0 $ 

25 

MPYAD 

WGNN,MJRRO,/MJRR $ 

26 

MATPRN 

MJRR, , , , $ 

27 

TRNSP 

MJRR/MJRRT $ 

28 

SOLVE 

WRD,MJRRT/MJRRDT/C,N,l $ 

29 

ADD 

MJRRDT,/MJRRPT/C,N, (100.,0.) $ 

30 

TRNSP 

MJRRPT/MJRRP $ 

31 

MATPRN 

MJRRP,,,,// $ 

32 

END ALTER 
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Explanation of DMAP Alter Statements 


I. Generation of Large Masses for SUPORT Points for Base Motion Analyses 

1. DMAP statements 2 through 14 are inserted after module RBMG4 
has run. For COSMIC Rigid Format 3 (April 1982 version) the 
Alter is at a^^= 69 and for MSC Rigid Format 25 (Dec. 1981, 
version 61B) a^ is 140. 

2-4. PARAM modules set the values of several parameters to be used. 

Parameter GRAV is the inverse of the regular WTMASS parameter. 

5. DIAGONAL extracts the diagonal of the rigid body mass matrix, 

WR ( equation (11)) and calls the result matrix data 

block MRD. 

6. Matrix data block WRD is MRD multiplied by parameter 
GRAV=1/WTMASS SO that WRD is generally in weight units if a 
PARAM WTMASS card is used in the bulk data. 

O 

7. Matrix data block MBASE is 10 times MRD. These give the 
large fictitious masses which are put or the r DOF for base 
motion problems. 

8-9. Module MERGE, using partitioning vector CP generated in module 
VEC, expands the large fictitious mass matrix, MBASE, up to 
full "a" set size. Output matrix MAAB has the large fictitious 
masses for the r DOF but in a matrix of "a" size. 

10. The large masses are printed out and can be seen to only exist 
at the r DOF. 

11-12. The large masses are added to the normal mass matrix MAA and 

the result equivalenced back to MAA so that when MAA is called 
for in later modules the sum of the normal MAA and the large 
masses MAAB will be used. 

13-14. The real eigenvalue analysis module, READ, uses the rigid body 
mass matrix, MR, when there are SUPORT DOF. Since the large 
masses were added at the r DOF, the sum of the MR (from the 
finite element model) and the large masses (in MBASE) is to be 
used for MR in module READ. 


II. Generation of Modal Participation Factors 

15. Following module READ and the output of the eigenvalues and 
eigenvalue analysis summary, the statements 16 through 23 are 
inserted for calculation of the modal participation factors. For 
COSMIC Rigid Format 3 (April 1982) the alter is at a =77 and for 
MSC Rigid Format 25 (Dec. 1981 version 613) a^ is 146. 
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16-17. Module PARTN extracts the eigenvectors for the il DOF and module 
GKAM extracts from this the elastic I ) eigenvectors. The 
result is called PHILN and is matrix used in equation 

(12). Matrix data block MGNN is the matrix of generalized 
masses, in equation (14). 

18, Convert MGNN to weight units if a PARAM WTMASS bu]k data card 
is present. 

19-20. Calculate r Mu) and call the result matrix data block 

ML3. Note that matrix data block DM is matrix Dji^-. Matrix 
data blocK WL3 is ML3 in weight units if PARAM WTMASS is 
present in cu.lk data. 

\r "■ ^ 


21-23. Calculate 


•ik- 


and print out C » the matrix of modal participation factors. 

III. Statements 24 through 31 calculate and print the modal masses. These 

cannot currently be calculated in COSMIC NASTRAN since the COSMIC version 
of module DIAGONAL does not have a feature which allows calculation of a 
matrix which contains the terms of some input matrix raised to some 
specified power. 

24. Module DIAGONAL calculates a matrix, MJRRO, whose terms are the 
squares of the terms in the modal participation factor matrix, 
FNR. 


0 * 


25-26. By multiplying MJRRO by the generalized mass (or weight) matrix 
the diagonal terms of the modal mass lattices are obtained. 
Matrix is so arranged that jts j row is the diagonal 

of the j modal mass matrix, , in equation (27) . Tl.is 

matrix is printed out. 

27-31. The modal mass terms are divided by the corresponding rigid 
body mass terms for the complete structure and multiplied by 
100. The resulting matrix, MJRRP, is printed out and gives the 
diagonals of the modal masses (in the rows of MJRRP) in percent 
of the rigid body mass. If there are no masses (other than the 
large fictitious masses) at the r DOF then the columns of MJRRP 
should add to 100 if all modes are used. 
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